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Abstract 

Limits on dark matter spin dependent elastic scattering cross section on protons 
derived from IceCube data are obtained for different dark matter annihilation chan¬ 
nels using micrOMEGAs . The uncertainty on the derived limits, estimated by using 
different neutrino spectra, can reach a factor two. For all dark matter annihilation 
channels except for quarks, the limits on the spin dependent cross section are more 
stringent than those obtained in direct detection experiments. The new functions 
that allow to derive those limits are described. 
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1 Introduction 

The strong astrophysical and cosmological evidence for dark matter (DM) motivates nu¬ 
merous direct and indirect searches for DM both in astroparticle experiments and at 
colliders. Exploiting eventual signals from different sources could shed light on the true 
nature of DM. Indirect searches for DM annihilation in the galaxy are actively pursued 
using either positrons or antiprotons [1], photons [2] and neutrinos [3,4]. The best limit 
on galactic neutrinos have been achieved by IceCube [4], however the limit on the DM 
annihilation cross section are still orders of magnitudes above the canonical cross section 
required to achieve the relic density assuming a standard cosmological scenario. Neutrino 
telescopes such as Super-Kamiokande [5], Baksan [6], Amanda and IceCube [7] can also 
observe neutrinos originating from annihilation of dark matter captured in the Sun. In 
that case, the neutrino flux is determined by the cross section for DM scattering on nuclei 
which drives the capture rate and is thus related to direct detection searches. Neutrino 
telescopes are sensitive to both DM nuclei spin dependent (SD) and spin independent 
(SI) interactions. However SI interactions are currently strongly constrained by direct 
detection experiments such as LUX [8] and XENON [9] whereas there is much more free¬ 
dom for the SD case [10]. Moreover the coherent enhancement of SI interactions only 
occurs for heavier nuclei which are subdominant in the Sun. The prospects to constrain 
the parameter space of DM models such as the minimal supersymmetric standard model 
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(MSSM) using the neutrino flux from DM capture in the Sun with IceCube was exam¬ 
ined in [11]. Furthermore limits on the spin dependent cross section based on IceCube79 
results have been compared to those obtained at the LHC from monojet searches within 
the framework of effective field theory for a choice of DM - quark effective operators [12], 
neutrino telescopes were shown to be more sensitive than colliders for heavy DM masses. 

In this article we present new features of micrOMEGAs [13-15] that allow to derive 
the bounds on spin dependent interactions of DM with protons using the data from 
IceCube22, and to obtain a likelihood that allows to combine these results with those 
from other searches when scanning over the parameter space of a generic DM model. We 
derive the bounds on spin dependent interactions with protons for any DM annihilation 
channel in standard model (SM) particles and compare these results with those obtained 
by the IceCube collaboration in two specific channels. We use the publicly available data 
of IceCube22, which gives a limit better than expected from statistics only because of a 
deficit in the number of observed events at small angles. For DM masses above 250GeV , 
the limits thus obtained are expected to be comparable to those that can be reached using 
the newer data from IceCube79 [16]. However for lighter DM, the lower energy threshold 
of IceCube and DeepCore leads to a higher sensitivity than IceCube22. In the case where 
DM annihilates into quarks, we compare our results with those derived from monojet 
searches at the LHC and interpreted within an effective field theory approach assuming 
the same spin dependent effective operator used for DM capture. We furthermore examine 
the impact of IceCube22 data on the parameter space of dark matter models, for this we 
use a simple Z-portal model as well as U(l) extensions of the MSSM. 

This article is organized as follows. Section 2 summarizes the equations used for 
computing the neutrino flux, section 3 describes the statistical analysis used to extract 
a limit. Section 4 presents the results for the upper limit on the spin dependent cross 
section for each individual channel. Section 5 shows the impact of IceCube results on the 
parameter space of two sample dark matter models. Our conclusions are summarized in 
section 6. The micrOMEGAs functions are described in the appendix as well as the method 
used to compare our results with DarkSUSY [17]. 


2 Neutrino Flux 


After being captured, DM particles concentrate in the center of the Sun and then annihi¬ 
late into Standard Model particles. These SM particles further decay producing neutrinos 
that can be observed at the Earth. The equation describing the evolution of the number 
of DM particles N x (assuming the DM is self-conjugate) reads 

N x = C x — A XX N X — EN X , (1) 


where A xx is the rate of depletion of DM, 


A 


xx ~ 


M xx 

Veff 


( 2 ) 


( av )xx the velocity averaged annihilation cross section of DM into SM particles and 
V e ff is the effective volume of DM in the Sun. E is the evaporation rate and C x is the 
capture rate of DM particles in the core of the Sun. It depends on the DM - nucleus 
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scattering cross section, as well as on the DM velocity distribution and local density [18]. 
The computation of the capture rate in micrOMEGAs is described in [14] and involves 
a contribution from spin independent and spin dependent interactions of DM on nuclei. 
The spin independent interactions add coherently, so the contribution of heavy nuclei 
is enhanced with a factor oc A 2 . However spin independent interactions are strongly 
constrained from direct detection experiments and therefore can safely be neglected as 
will be demonstrated in section 4. For spin dependent interactions there is no coherence 
effect and hydrogen gives the largest contribution, neutrino telescopes are therefore mostly 
sensitive to spin dependent interactions on protons. 

Neglecting the evaporation rate, the solution of Eq. (1) is simply 

r„ = \a xx n 2 x = tan! qycSv) (3) 

where t = 4.57 x 10 9 yrs is the age of the Sun. When the capture rate and the annihilation 
cross section are sufficiently large (y/C x A xx t 1) equilibrium is reached and the anni¬ 
hilation rate only depends on the capture rate, T xx = C x j 2 and not on DM annihilation 
properties. 

The total neutrino (antineutrino) spectrum at the Earth is proportional to the capture 
rate but also features a dependence on the DM annihilation channel 


dj(j)iy 

dE v 


^d 2 J xx f 




d_Nf_ 
ff dE v 


(4) 


where d Q — 1.5 x 10 8 km is the distance to the Sun, Br^j the branching fraction into each 
particle/antiparticle final state ff and Nf are the neutrino spectra resulting from those 
annihilations. The neutrino spectra originating from different annihilation channels into 
SM particles and taking into account oscillations and Sun medium effects were computed 
in WimpSira [19], DMz/ [20] and PPPC4DMz/ [21] 1 . We use the set of tables provided by 
these groups. 

The distribution of the number of events in the detector, N s , will depend on the flux 
of muon neutrinos. After an exposure time t e the number of events is given by 


dN s 

— £ = t e 
dE 



Ay{E) + 



(5) 


where A V (E) and A P (E) are the effective area of the detector for neutrinos and anti¬ 
neutrinos respectively. 


3 Limits on the number of signal events 

To obtain the limits on the neutrino flux that could originate from DM capture in the 
Sun, one first needs to distinguish the signal and background events. The neutrino fluxes 
depend on the angle ip between the incoming neutrino and the direction of the Sun. 
Background events are expected to be distributed over any angle ip whereas all signal 

■'■Note that neutrino signals from the Sun and corresponding limits on spin dependent interaction cross 
sections were also investigated in [22] by focusing on the effect of internal bremsstrahlung and in [23] by 
paying attention on consequences of final state radiation of electroweak gauge boson. 
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events should occur at ip = 0. However the IceCube experiment has a finite angular 
resolution, the angular distribution of signal events is therefore 


dN 


cos ip— 1 

= Ce-^. 


~~ , ( 6 ) 

dcostp 

where C is a normalization constant and a is the mean angle (in radians) [24] . The mean 
and median angle, p me d (in degrees) are simply related through 

[Pined 


7r 


cr = 


V^ 180 ' 


(7) 


The energy dependence of the median angle is provided on the IceCube official website - ’. 
For example, for a DM of mass 1000 GeV which annihilates into W + W~, p me d = 2.9° [25]. 

After estimating the background using neutrino events registered at large angles, an 
upper limit on the number of signal events can be obtained from the number of events 
registered at small angle, this method allows to get rid of systematic uncertainties. To 
obtain an upper limit on the number of signal events, IceCube22 uses the Feldman-Cousins 
method [26]. In IceCube22, 13 events were registered in the angular bin ip < 3° while 18 
background events are expected. From the number of events in that angular bin we have 
checked that the upper limits on the total number of signal events, /x 90 , agree with the 
ones presented in [25] and displayed in the second column of Table 1 for the case of DM 
annihilating into W + W~. 


rn x 

fi-90 [2 j] 

micrOMEGAs 

PPPC4DM^ 

micrOMEGAs 

WimpSim 

DarkSUSY 

IceCube22 

. .micrO 
^90 

250 

7.5 

1.37 10 -40 

9.51 10~ 41 

8.9 10 -41 

2.8 10- 40 

9.07 

500 

6.8 

1.48 10“ 40 

1.04 10“ 40 

9.7 10“ 41 

3.0 10“ 40 

8.98 

1000 

6.8 

4.38 10" 40 

2.67 10~ 4 ° 

2.5 10" 4 ° 

8.7 10" 40 

8.95 

3000 

6.4 

3.00 10" 39 

2.06 10" 39 

2.1 10" 39 

9.9 10" 39 

8.65 

5000 

6.8 

7.84 10" 39 

6.00 10" 39 

5.8 10" 39 

3.6 10" 38 

8.67 


Table 1: 90%CL upper limit on the number of signal events, /i 90 , for DM of mass m x in 
[GeV] annihilating into W + W~ and the corresponding limits in [cm 2 ] on spin dependent 
DM - proton cross section obtained by micrOMEGAs with PPPC4DMz/ or WirapSim , by 
DarkSUSY , and by IceCube [25]. The last column contains the 90%CL upper limit on the 
number of signal events /rgo evaluated using the likelihood function of micrOMEGAs . 


An alternative approach to derive the upper limit on the number of events based on 
Bayesian statistics has been implemented in micrOMEGAs . The likelihood function uses 
both information about the angular distribution and the number of digital optical modules 
(DOM) activated for each event. Note that the number of active DOMs is correlated 
with neutrino energy. The method makes use of the experimental data provided by the 
IceCube22 collaboration and in the DarkSUSY-5.1.2 distribution [24] 2 3 , this includes, 

2 https://icecube. wisc.edu/science/data/ic22-solar-wimp 

3 shared/DarkSusy/IC_data 
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• A U (E) - the effective area for neutrino detection; 

• Ap(E) - the effective area for anti-neutrino detection; 

• <j(E) - the mean value of the neutrino angle resolution; 

• chk(N c han) - the probability to activate N c h an DOMs [24], This is given for 20 energy 
intervals in the range 39 GeV < E u < 4 10 5 GeV; 

• , djVbg - the distribution of background events as a function of the angle between the 
incoming (anti)neutrino and the Sun; 

• chb g (N c han) - the probability to activate N c h an DOMs for background event; 

• t e the exposure time t e = 104/365 [Year], 


For background events, the angular and energy distribution are not expected to be 
correlated. Then the probability of observing a background event with an angle ip that 
activates N chan DOMs in the detector is 


Pfr-y (cOS ip. N c h an ) 


dN bg 
d cos 


ch bg (N c 


chan ) 


( 8 ) 


The probability of a signal event for given neutrino </F and anti-neutrino fluxes 
is defined as 


20 


where 


^ 77-2* ' 

Ps (cos ip, N chan ) = 2^ — e 
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2 
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n k Je dE 
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( 10 ) 

( 11 ) 


The IceCube collaboration also provides the full list of events registered during the 
exposure time t e , altogether 6946 events. Most of these events occur at large angle with 
respect to the direction of the Sun. To construct a likelihood function we avoid the regions 
with a large atmospheric background, and we will consider only events at angles less than 
ip cut . The number of expected events for both background (N bg ) and signal (N s ) is then 


COS (pcut 


N bg = 5Z / ^bg(cOSip,N chan )dcOSip 

^chan Q 

COS (flout 

-V = y" j V,(cosip,N Am )dcosip. 


( 12 ) 

(13) 


The likelihood function reads 

£( a , F cut ) 


(N? +aN c ) Ntot 

K^bg ' ucly s) -jV£-aiVf 


Ntot 


N 


-e b s 


tot * 


IP(cOS (pi , N c han,ii ^0 


(14) 
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where N tot is the total number of detected events for p < tp cut . The coefficient a varies 
in the range [0, oo] and 

F(cOSCp,N chan ,a) ' ~— i^bgi^COS (p, N c h an ^ T ttlP s (COS <p, Nchan )) ■ (15) 

N bg + aN s 

To obtain the limit on the number of signal events we use the Bayesian approach with a 
flat prior distribution for the number of signal events. Then the confidence level or rather 
credibility level, is given by [27] 


CL 


1 


OO 


I C(a)da/ 
o 


/ 


C(a)da, 


(16) 


Here CL represents the probability that the number of signal events is less than N r s . 



^cut 


Figure 1: 90%CL limit on the total number of signal events as a function of tp cut for 
a DM mass of 1 TeV in the W + W~ channel obtained with micrOMEGAs (black) and 
DarkSUSY (dash/red). 


Fig. 1 shows the dependence of the 90%CL upper limit on the number of signal events 
on <p cut . Clearly the limit does not depend on p cut for p cut > 5°, in our calculations we fix 
'Lent = 10°. For comparison the dependence on p cut of the upper limit on the number of 
signal events obtained with the default option in DarkSUSY is also displayed. Note that the 
two methods give similar results for ip cut = 3.7°, this angle also corresponds to the maximal 
value for N^/y/N£ ckg . The 90%CL limits on the number of signal events calculated by 
micrOMEGAs are presented in the last column of Table 1. Despite the different statistical 
methods used, these numbers are in rough agreement with those extracted using only one 
angular bin as done by IceCube22 [25]. 

Note that to take into account the 20% uncertainty on the effective area, micrOMEGAs 
divides A„(E), Ap(E) by 1.2 when calculating the CL. Therefore we effectively rescale by 
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20% the limit extracted on the SD cross section and we thus obtain a more conservative 
limit. Note that this rescaling has no impact on the upper limit on the number of signal 
events. 


4 Extracting limits on DM — proton spin dependent 
cross sections 

For a given number of signal events, one can extract an upper limit on the DM - nuclei 
cross section, this limit will depend on the fluxes of muon neutrinos/antineutrinos that 
reach the detector, Eq. (4). This spectrum will depend primarily on the DM annihilation 
channel. When neutrinos propagate through the Sun, the spectra will be distorted by 
neutrino attenuation in the Sun and neutrino oscillation, furthermore other effects like 
electroweak bremstrahlung can also play a role. The latter is particularly important for the 
electron and muon channels, otherwise these channels would not lead to neutrinos because 
the electron is stable and the muon is quickly absorbed in the Sun. Two sets of neutrino 
spectra are publicly available, WirapSim [19] and PPPC4DMz/ [21], Both spectra are 
implemented in micrOMEGAs and are used to compute the neutrino flux. We have checked 
the validity of our implementation by comparing with the spectra presented in [19, 21] 
and directly with the WimpSim tables implemented in DarkSUSY-5.1.2 [17]. Furthermore 
we have also compared our current implementation with the previous version of neutrino 
spectra DMz/ [20] implemented in micrOMEGAs [14]. 
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Figure 2: Comparison of neutrino spectra from PPPC4DMz/ (dashed/blue) and 
WimpSim (full/red) for DM annihilation into bb (left) and W + W~ (right). The DM mass 
is set to 500 GeV for bb and 1 TeV for W + W~. The spectra from the DMz/ tables [20] 
are also displayed for comparison (dotted, light green). 


For some channels the WimpSim and PPPC4DMz/ spectra can be quite different, see 
Fig. 2 which compares the total neutrino/anti-neutrino spectra for two different channels. 
This difference will impact the limits derived on the DM - nucleon cross sections. The 


7 




reason for these differences is not completely clear despite the different approaches used. 
WimpSim uses Pythia [28] to get the initial neutrino spectra and a Monte Carlo approach 
for neutrino propagation while PPPC4DMz/ uses GEANT [29] for the initial spectra and 
solves the equation for the density matrix. The discrepancy in the bb channel could be 
due to differences in Pythia and GEANT since the previous version of the neutrino spec¬ 
trum DMz/ [20], which relied on Pythia is in perfect agreement with WimpSim . We have 
also noted that the different treatment of the propagation in the Sun leads to different 
neutrino spectra even in the case where DM annihilates directly into neutrino pairs, only 
the WimpSim spectrum featuring a peak at the maximum energy. In such case it is not sur¬ 
prising that for all channels with an initial hard neutrino spectrum, for example W + W~, 
WimpSim leads to a larger and harder flux, see the right panel of Fig. 2. Another differ¬ 
ence between the two codes is that bremstrahlung is included in PPPC4DMi/ , however 
bremstrahlung is expected to be important mainly for light leptons. 

The neutrino spectra produced by polarized vector bosons and leptons are also pro¬ 
vided in PPPC4DM^ . The polarisation of vector bosons affects the spectrum of high en¬ 
ergy neutrinos at the source, after propagation the effect is of the order of 10%. Although 
mi crOMEGAs has an option to automatically compute the polarization of vector bosons 
in DM annihilation processes, the neutrino signal in a given model is always computed 
assuming unpolarized vector bosons since the effect is not large. The polarized tables 
can be used only in a model independent approach when one inputs a cross section and 
chooses explicitly the annihilation channel. On the other hand micrOMEGAs does not keep 
track of the polarization of leptons. The polarization effect can be much more important 
for annihilation channels into light leptons, with in particular an important contribution 
from left-handed light leptons (e, /i) due to bremstrahlung. However the neutrino flux is 
expected to be very small for these channels as will be demonstrated in the next section. 
Furthermore the neutrino spectra only depend weakly on the r polarization. 

4.1 Comparison with DarkSUSY 

The value of the DM - proton cross section corresponding to the upper limit on the 
number of signal events follows from Eqs. (3), (4) and (5). The results computed with 
micrOMEGAs and using two different choices of neutrino spectra [19], [21] are also presented 
in Table 1. These values are obtained starting with the upper limit on the number of 
events given in the first column and obtained in [25]. The results obtained with the 
WimpSim spectra are compared with those obtained with DarkSUSY [17] and are found to 
be in excellent agreement - for a description of the method used see appendix C. The limits 
presented by the IceCube collaboration on the spin dependent cross section are higher 
by roughly a factor of three [25] and are consistent with the ones obtained using the 
DarkSUSY option to compute the 90%CL using (p cut = 8°. Note that the limits obtained 
with DarkSUSY strongly depend on the value chosen for the angle. Since the signal is 
concentrated at small angle, the limits can be improved with a smaller angular cut. To first 
approximation, the optimal cut corresponds to the maximum of N S i gna i (4> cu t) / yf A T b g ( (j) cu t ) • 
For a 1 TeV DM annihilating into W + W~, the maximum is reached at (p cut = 3.7°. For 
this cut, the 90% exclusion limit on the number of signal events obtained with DarkSUSY is 
8.5, in good agreement with the micrOMEGAs result shown in the last column of Table 1. 

In Table 2 we compare the 90% confidence limit on the DM-proton SD cross section ob¬ 
tained with DarkSUSY and with micrOMEGAs for the case of DM annihilation into W + W~. 



For this we use the default “Number likelihood” flag of DarkSUSY , For each DM mass 
we find the angular cut which minimizes the ratio N signal (<j ) cut )/ y/N bg (<f) cut ) where N signa i 
is the predicted number of signal events and iV& s is the number of background events 
estimated by DarkSUSY . Note that 4> cu t does not depend on the DM-proton cross section 
nor on data. These results are compared with the ones obtained with mi crOMEGAs using 
only angular distribution. The value obtained with micrOMEGAs is larger because we allow 
an additional 20% uncertainty to take into account the systematic error on the effective 
area of IceCube22. We also compare with the results obtained with micrOMEGAs adding 
the information from the energy distribution as described in the previous section. Here 
limits are improved noticeably for DM masses below 1 TeV, while at higher masses the 
two methods give similar results, this is because the energy distribution of signal and 
background events are similar. The same comparison is performed for the bb channel and 
presented in Table 3. 


m x 

(GeV) 

IceCube22 

[25] 

DarkSUSY 

with 0 cut 

micrOMEGAs 

0 only 

micrOMEGAs 

0 N c han 

IceCube79 

[16] 

100 


4.19 10 _39 (5.6°) 

4.65 10- 39 

3.25 10- 39 

2.68 10~ 40 

250 

2.8 10~ 40 

1.33 10 _40 (4.3°) 

1.56 10- 40 

1.38 10- 40 

1.34 10~ 40 

500 

3.0 10 -40 

1.67 10 _40 (3.9°) 

1.71 10“ 4 ° 

1.64 10- 40 

1.57 10~ 4 ° 

1000 

8.7 10 -40 

3.12 10 _40 (3.7°) 

4.13 10- 40 

4.21 10- 40 

4.48 10- 40 

3000 

9.9 10" 39 

2.73 10 _39 (3.7°) 

3.32 10- 39 

3.32 10" 39 

5.02 10" 39 

5000 

3.6 10" 38 

7.24 10 39 (3.7°) 

9.09 IQ" 39 

9.09 10" 39 

1.59 10~ 38 


Table 2: Comparison between limits on a xg (cm 2 ) for xx W + W~ obtained by IceCube, 
DarSUSY for the value of (f) cut given in parenthesis and micrOMEGAs using only angular 
information (0 only) and also energy distributions (0 + N chan ). All results are obtained 
using the WimpSim tables. 


Finally, we comment on more recent results from the lceCube79 collaboration [16]. 
The lower energy threshold achievable with the use of DeepCore leads to a significant 
increase in sensitivity for DM masses below ~ 250 GeV. At high masses the limits that 
are obtained by the collaboration in the W + W~ and bb channel are roughly a factor 2 
better than those of IceCube22 [25], see Tables 2 and 3. However we have already shown 
that the likelihood function implemented in micrOMEGAs leads also to better limits than 
those found in [25]. In part it is due to a statistical anomaly associated with the small 
number of detected events at small angles (13 events observed and 18 events expected). 

Note that using the PPPC4DMz/ spectra leads to a more conservative limit in the WW 
channel, this is expected considering the difference in neutrino spectra mentioned above. 
This is not the case in all channels as will be discussed in the next section. 
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m x 

(GeV) 

IceCube22 

[25] 

DarkSUSY 

with <p cut 

micrOMEGAs 

</> only 

micrOMEGAs 

4> ~\~ Nchan 

IceCube79 

[16] 

100 


1.63 10 _36 (6.2°) 

3.74 10“ 36 

2.59 10~ 36 

1.47 10~ 38 

250 


2.94 10~ 38 (5.0°) 

3.23 10“ 38 

2.40 10~ 38 

5.90 10~ 39 

500 

2.6 10“ 38 

1.13 10 _38 (4.5°) 

1.13 10“ 38 

9.46 10- 39 

7.56 10~ 39 

1000 

2.2 10“ 38 

1.02 10 _38 (4.2°) 

1.16 10“ 38 

1.05 10~ 38 

1.00 10~ 38 

3000 

7.2 10" 38 

2.69 10 _38 (4.0°) 

3.55 10~ 38 

3.35 10" 38 

3.16 10" 38 

5000 

1.5 10" 37 

4.75 10~ 38 (4.1°) 

7.43 10" 38 

7.28 10" 38 

7.29 10" 38 


Table 3: Same as Table 2 for XX bb. 


4.2 Model independent limits on DM — proton cross section 

Using the likelihood method described in section 3, we can extract the limits on spin 
dependent interactions of DM on protons for each channel of DM annihilation into SM 
particles for both the PPPC4DMi/ and WimpSim spectra. The results are displayed in 
Fig. 3 and Fig. 4. 

The strongest limits are obtained when DM annihilates into neutrinos and in particular 
into Up or v T . We find a SD > 9.7 x 10~ 42 pb for m x = 250 GeV using WimpSim . At masses 
above the TeV scale the v T channel becomes much more sensitive than the channel. 
This is due to the fact that v T can be regenerated, a charged current interaction of v T 
with nuclei generates a r which in turns decay rapidly into v r with a high probability. 
On the other hand the /i generated from by charged current interactions is stopped in 
the Sun medium before emitting a neutrino while the electron generated from v e is stable. 
Note that the different treatment of propagation in WimpSim and PPPC4DMzz leads to 
roughly a factor 2 in the 90%CL limit on , with the PPPC4DMi/ leading to weaker 
limits because the neutrino spectrum is softer. For all masses above 60GeV the neutrino 
telescope limits are more stringent than the ones from direct detection. 

At high masses the second most sensitive leptonic channel is the r channel, furthermore 
for all masses above 100 GeV sensitivity is better than direct detection. Note that DM 
annihilating into fi or electron pairs can also lead to a neutrino flux, even though muons 
are quickly absorbed in the Sun and electrons are stable. This is due to electroweak 
bremstrahlung [21]. However the limits obtained are much weaker than for the t + t~ 
channel. The effect of including lepton chirality is illustrated for the /i + /r~ channel and 
basically leads to an improvement on the limit on SD cross section by a factor 2, this 
channel gives slightly more stringent limits than PICO in the range 600 - 1000 GeV. 

The results for the W/Z and Higgs channels are displayed in Fig. 4 (left panel). The 
most sensitive channels are respectively ZZ and W + W~ for DM masses below the TeV 
scale, while above 2 TeV the channel HH becomes slightly more sensitive. In fact at 
masses of 5 TeV the vector boson channels are within a factor 3 above the best leptonic 
channels. Using the PPPC4DMi/ tables we have extracted independently limits from the 
Wt and Wl channels, the results are within 10% of those for the unpolarized case. We have 
also computed the limits on a SD assuming that DM annihilates solely into gluon or photon 
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m x (GeV) m x (GeV) 


Figure 3: Upper limit on for DM annihilation in z/yV i = e, /j,,t (left panel) or t + t~ 
(right panel) using the PPPC4DMz/ (full) or WimpSim (dashed) spectra. The best limit 
on from the direct detection experiment PICO [10] is shown for comparison. The 
right panel also shows the limits from the light lepton channels using PPPC4DMz/ (full). 
The limits for polarized muons (dashed) are also displayed. 


pairs, however they are weak. We find that for m x = 800 GeV, a SD > 6.6 x 10 -38 pb in 
the gluon channel and a SD > 1.9 x ICC 38 pb in the photon channel. When DM annihilates 
into quarks the best limits are obtained for heavy quarks, in particular from the tt channel 
which lies about a factor of 3 above that of the ZZ case. In all those limits there is clearly a 
large uncertainty coming from the computation of the neutrino spectra as discussed above. 
Note that for light quarks it is the PPPC4DMz/ spectra that leads to the most stringent 
limits, as expected since in this case the neutrino flux is larger, see Fig. 2 for the bb 
channel. In fact for DM above 2 TeV, the bb channel has a comparable sensitivity (within 
a factor two) to the top or boson channels. For both the third generation quarks the 
upper limit on is better than the one obtained in direct detection for masses above 
a few hundred GeV’s. For light DM masses, these limits are however weaker than those 
obtained from the ATLAS and CMS collaborations in searches for DM in monojet events 
with 8 TeV data at the LHC [30,31]. The collider limits presented in Fig. 4 are based on 
an effective field theory (EFT) approach and assume that DM couples to all quarks via 
the operator l/M^x'y^^XQlfi'y 5 ■ I n this approach the production of DM at colliders can 
be directly related to the direct detection cross sections provided the scale of new physics 
inducing this operator is large enough for the EFT to be valid. 

We have also extracted the limits on a SI assuming DM annihilates completely in v T 
pairs, here we assume equal cross section for scattering on protons or neutrons. We find 
that even with this favorable channel, neutrino telescopes cannot compete with direct 
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m x (GeV) m x (GeV) 


Figure 4: Upper limit on cr^f for DM annihilation in W + W~, ZZ , HH, 77 (left panel) or 
bb,cc,tt,qq (right panel) using the PPPC4DMz/ (full) or WimpSim (dashed) spectra. The 
best limit from the DM direct detection PICO [10] is also displayed as well as the limits 
from mono jet analysis in CMS [30] and ATLAS [31] 


detection experiments, the best limit being a factor 2-10 above that of LUX [8], see 
Fig. 5. 

5 Constraints on DM models 

In a typical DM model, DM annihilation can proceed through different channels. We 
investigate the constraints that can be derived from IceCube data within two sample 
DM models, a Z-portal model [32] and U(l) extensions of the MSSM with a neutralino 
DM [33,34]. Note that for this analysis, we include both SD and SI interactions and that 
we use the WimpSim spectra. 

5.1 Z-portal dark matter 

This model is a simple extension of the SM with a fermion dark matter candidate 
which couples to the Z with either vector ( V x ) or axial-vector couplings (A x ), [32], 
Such interaction can be obtained after symmetry breaking from the effective operator 
ig/v 2 D(xi 11 (V x + A x 75 )%) where H is a scalar doublet. The relevant DM interac¬ 
tions are described by the Lagrangian 

(IT) 
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m x (GeV) 


Figure 5: Upper limit on cr^ for DM annihilation into v T v T using the PPPC4DMz/ (full) 
or WimpSim (dashed) spectra compared with the best limit from DM direct detection in 
LUX [8], 


In addition to the DM mass the model contains only two free parameters, we choose 
A x and the ratio a = A x /V x . The SD interaction which proceeds through Z exchange 
depends only on the axial-vector coupling (A x ) while the SI interaction depends on the 
vector coupling (I%). When kinematically accessible, DM annihilation is mainly into W 
pairs, which depends on V x , or into final states involving ZZ^ZH and heavy fermion 
pairs which mostly depend on A x . Annihilation into leptons and light quarks is generally 
suppressed. Note that terms involving the Higgs in Eq. (17) are necessary to insure 
gauge invariance and have to be include in the relic density computation. To assess the 
constraints originating from IceCube22, we have performed a random scan over the free 
parameters of the model in the range 

50GeV < m x < 5000GeV 10“ 5 < \A X \ < 1.0 0.1 < |a| < 10 5 . (18) 

We have imposed the relic density constraint from PLANCK, allowing a 10% uncer¬ 
tainty 0.9 < O/z 2 /0.1199 <1.1 [35]. For the allowed points, we then determined those 
that are excluded by LUX and/or by IceCube, see Fig. 6 . The former excludes points 
where V x is large while IceCube constrains points where A x is large. The DM mass is 
also relevant, for tudm < m w annihilation into light quarks dominates and this region 
is thus poorly constrained by IceCube. In the region where mw < »»dm < mtop we find 
the best IceCube exclusion since annihilation into gauge bosons leads to strong exclusion. 
Heavier DM annihilates rather in tt or ZH, both have weaker exclusion power than the 
W + W~ channel. For these masses it is generally LUX that excludes the points. After 
applying the PLANCK constraints, for a given dark matter mass, most points in Fig. 6 
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m x (GeV) 


Figure 6 : Constraints on the parameter space of the Z-portal model in the a^—m x plane. 
The color code indicates points that are constrained by IceCube22 and by LUX (green) 
by IceCube22 only (black) and by LUX only (red). Blue points satisfy all constraints and 
have a relic density in agreement with PLANCK. 


are concentrated in a narrow band for the SD cross section. This is somewhat an artifact 
of the scan, allowing very small values of a (or a dominant V x coupling) would populate 
the region below the band. However, a large value for V x also leads to a large spin inde¬ 
pendent cross section, thus all those points are constrained by the LUX result even if in 
some case they can be also constrained by IceCube. In summary after imposing the upper 
limit on SI interactions, IceCube can further constrain DM in the narrow mass range 80 
- 175 GeV. 

5.2 UMSSM 

U(l) extensions of the MSSM (UMSSM) are characterized by additional vector and scalar 
superfields. The model therefore features a new gauge boson ( Z') and a singlet scalar, 
each with their supersymmetric partner, furthermore the model contains right-handed 
(RH) neutrinos. The lightest neutral supersymmetric particle can be the DM candidate, 
either a nentralino or the partner of one of the RH neutrino. Moreover new tree-level 
contributions to the Higgs mass make it possible to achieve a Higgs mass at 125 GeV 
without large corrections from the stop sector. We use the results of a random scan 
performed in [34] where constraints from flavor observables, the Higgs mass and couplings, 
and LHC searches for the Z' and for supersymmetric particles were considered. Moreover 
an upper bound on the relic density of dark matter was imposed [35] as well as an upper 
bound from direct spin independent dark matter searches from LUX [8]. We then compute 
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the CL for IceCube exclusion for each point of parameter space. In computing the capture 
rate we took into account the fact that the lightest neutralino Xi could constitute only a 
fraction of the dark matter using a rescaling factor £ defined by 


0 ^ 68 for h 2 < 0.1168, 

1 for VL^oh 2 € [0.1168,0.1208], 


(19) 


where the 2cr deviation from the central value measured by PLANCK is f Vi 2 = 0.1168. 



Figure 7: Allowed parameter space of the UMSSM model in the £cr| P_ p — m^o plane. 


Our results for all points satisfying the above mentioned constraints are shown in 
Fig. 7 in the spin dependent cross section vs neutralino mass plane. We did not find any 
point falling in the 90%CL exclusion by IceCube22, although we found such points when 
we ignored the LUX constraint. Moreover when assuming some regeneration mechanism 
for DM, which implies ignoring the rescaling of the DM local density, a few points were 
excluded by IceCube. Since our limits are too conservative at low masses we expect 
that IceCube79 data would exclude some of these points. Indeed a prospective study 
had shown that in the MSSM several points were excludable by IceCube+DeepCore [11]. 
Moreover after applying the upper limit on the SI direct detection cross section such 
points are concentrated at masses below 250 GeV. Although the two studies are not 
directly comparable since they are based on different supersymmetric models and we 
include most recent constraints form the LHC on the Higgs and supersymmetric sectors, 
it is encouraging that they qualitatively agree. 

6 Conclusion 

We have shown that upper limits on spin dependent interactions of DM with protons 
obtained from IceCube22 data are stronger than those from direct detection in a variety 
of channels and for a large range of masses above 100 GeV. The best limits are obtained 
for DM annihilation directly into neutrinos, although the t + t~ , W + W~, ZZ channels also 
display good sensitivity especially for masses of a few hundred GeVs . Uncertainties from 
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the propagation and the oscillation of neutrinos in the Sun estimated by comparing the 
neutrino spectra computed by two different methods induce up to a factor two difference 
on the extracted limits on the DM proton spin dependent cross section. Other source 
of uncertainties affecting the limits derived from DM capture in the Sun arise from the 
DM local density, the DM velocity distribution, the escape velocity in the Sun, and the 
presence of a dark disc. These uncertainties were estimated to impact the capture rate 
by at most 50% [36]. Note however that direct detection limits are also sensitive to the 
uncertainties on the DM local density and velocity distribution. 

A direct comparison of the LHC and neutrino telescope limits on the spin dependent 
cross section is straightforward in an EFT approach when we assume that DM couples 
to quarks only. In this case we found that the collider limits were better for light quarks 
and for light DM. 

We have also shown that the IceCube22 results constrain some area of the parameter 
space of typical dark matter models, although the potential is much more limited when 
the current constraints on the DM spin independent scattering cross section on nucleons 
are taken into account first. 

The results presented here and implemented in micrOMEGAs_4.2 are based on the 
publicly available data from IceCube22, the micrOMEGAs facilities to compute limits on 
DM proton cross section from neutrino telescopes will be extended when more recent 
IceCube data becomes publicly available. 
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Appendices 

A micrOMEGAs functions for neutrino telescopes 

Details on all micrOMEGAs functions can be found in the manual contained in the man 
directory of each version. Here we present the functions to compute the (anti-)neutrino 
fluxes as well as the likelihood and the CL corresponding to the data of Icecube22 for 
DM nucleon cross sections defined in a model-independent way or within the context of 
a specific model. 

We have introduced a new global parameter WIMPS IM which allows to choose the 
neutrino spectra. The default value WIMPSIM=0 corresponds to the PPPC4DMz/ [21] 
spectra while WIMPSIM=1 corresponds to the WimpSim [19] spectra and WIMPSIM=-1 to 
the DMi/ [20] spectra implemented in previous micrOMEGAs versions. The parameter 
forSun=l is used to compute the signature of DM annihilation in the Sun and Mdm des- 
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ignates the mass of the DM particle. 

• basicNuSpectra(forSun,Mdm,pdg, pol, nu_spect, nuB_spect) 

calculates the and P /; spectra corresponding to one DM annihilation into a particle- 
antiparticle pair with PDG code pdg. Mdra is the DM mass. Note that this routine depends 
implicitly on the global parameter WIMPSIM. The parameter pol selects the spectra for po¬ 
larized particles available in PPPC4DM^ . pol=-l(l) corresponds to longitudinal (trans¬ 
verse) polarisation of vector bosons or to left-handed (right-handed) fermions, pol=0 is 
used for unpolarized spectra. When polarized spectra are not available, the unpolarized 
ones are generated irrespective of the value of pol. The resulting spectrum is stored in 
arrays nu_spect and nuELspect. 

• captureAux(fv,forSun,Mdm,cslp,csln,csDp,csDn) 

calculates the number of DM particles captured per second assuming the cross sections 
for spin independent and spin dependent interactions with protons and neutrons cslp, 
csln, csDp, csDn are given as input parameters (in [pb]). A negative value for one 
of the cross sections is interpreted as a destructive interference between the proton and 
neutron amplitudes. The first argument is the DM velocity distribution, usually we 
substitute fv = Maxwell. The result depends implicitly on the global parameters rhoDM 
representing the local DM density, by default rhoDM = 0.3GeV/cm 3 . 

• IC22nuAr(E) 

effective area in [km 2 ] as a function of the neutrino energy, A v (E). 

• IC22nuBarAr(E) 

effective area in [km 2 ] as a function of the anti-neutrino energy, A p (E). 

• spectrMult(Spec, func) 

allows to multiply the spectrum Spec by any energy dependent function func. 

• spectrlnt(El,E2,Spec) 

integrates a spectrum/Flux, Spec from E\ to E 2 . 

• IC22BGdCos(cs) 

angular distribution of the number of background events as a function of cos<^, ■ 

• IC22sigraa(E) 

neutrino angular resolution in radians as a function of energy, Eq. ( 6 ). 

• exLevIC22( nu_flux, nuB_flux,&B) 

calculates the credibility level for number of signal events generated by given and 
fluxes, Eq. (16). The fluxes are assumed to be in [GeV km 2 Year] -1 . This function uses 
the IC22BGdCos(cs) and IC22sigma(E) angular distribution for background and signal as 
well as the event hies distributed by IceCube22 with p < p cut = 8°. The returned param¬ 
eter B corresponds to the Bayes factor, £(1, yw)/£(0, p cu t), with C defined in Eq. (14). 

• fluxFactorIC22(exLev, nu,nuBar) 

For given neutrino, nu, and anti-neutrino, nuBar, fluxes, this function returns the factor 
that should be applied to the fluxes to obtain a given exclusion level exLev in exLevIC22. 
This is used to obtain limits on the SD cross section for given annihilation channel. 

• IC22events(nu,nuB, phi, &Nsig,&Nbg,&Nobs) 

For given neutrino and anti-neutrino fluxes this routine calculates the expected number of 
signal Nsig, background Nbg and observed Nobs events, phi is the angular cut defined in 
degree units. The last parameter has to be an int type variable. Note that in calculating 
exclusion limits on cross sections, a factor 1.2 to take into account the uncertainty on the 
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effective area of the detector is included. 

To calculate the neutrino fluxes within some theoretical DM model one can use instead 
of basicNuSpectra and captureAux, the function 

• neutrinoFlux(fv,forSun,nu, nu_bar) 

which calculates the muon neutrino/anti-neutrino fluxes near the surface of the Earth. 
This function a) calculates the capture, annihilation, and evaporation rates in Eq. (1) and 
b) solves it numerically; c) calculates all branchings of DM annihilation and substitute 
them in Eq. (4) to get the fluxes. Here fv is the DM velocity distribution normalized 
such that / 0 °° vf(v)dv = 1. The units are km/s for v and s 1 /km 2 for f(v). For example 
one can use the Maxwell function. The calculated fluxes are stored in the arrays nu and 
nu_bar, the units used are [GeV km 2 Year]^ 1 . 

In micrOMEGAs the neutrino spectra and fluxes are tabulated in arrays of double nu¬ 
merical type with NZ = 250 elements. The Erst (zeroth in C) element of this array 
contains the maximal energy of the spectrum. It usually corresponds to the DM mass, 
see [14], The functions to interpolate and display spectra include 

• SpectdNdE(E,Flux) 

interpolates the tabulated spectra/flux and returns energy distribution. 

• displaySpectra(title, Emin, Emax, N, nul,labl,...) 

displays several spectra. Here title contains some text, Emin,Emax are the lower and 
upper limits, and N is the number of spectra to display. Each spectrum is defined with 
two arguments, nul designates the spectrum array and labl contains some text to label 
the spectrum. 


B Examples 

An example of the code that computes the 90% exclusion limits for DM-proton spin 
dependent cross sections for W annihilation channels and reproduces data of Tables 1, 2 
and 3 reads : 

double nu[NZ],nuB[NZ]; // arrays for neutrino and antineutrino spectra 
double csSDp [6]={1.E-40,9.51E-41,1.04E-40, 2.67E-40,2.06e-39,6.0e-39>; 

// cross section TAB_l,4-th column 
double dSun=150E6; // distance to Sun [km] 

double yrs=31556925.2; // year [sec] 

double Mdm[6]={100,250,500,1000,3000,5000}; // DM masses 
int forSun=l; 

double exLev=0.9; // 907, exclusion confidence level 
double Nsig; // for number of signal events 

WIMPSIM=1; 

printf("W channel:\n"); 

for(int k=0;k<6;k++) 

{ 

double Crate=captureAux(Maxwell,forSun,Mdm[k],0,0,csSDp[k] *1E36,0); //DM capture rate[s] 
double gamma=Crate*yrs/(4*M_PI*dSun*dSun)/2; 

double f_90; // improving factor for 907 confidence level exclusion 
printfC'Dark matter mass is 7.2E\n" ,Mdm [k] ) ; 

basicNuSpectra(forSun,Mdm[k],24/*W*/, 0, nu, nuB); // neutrino spectra 
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// neutrino fluxes 


for(int i=l;i<NZ;i++) { nu [i]*=gamma; nuB [i]*=gamma;} 
if (k) 

{ IC22events(nu,nuB, 10, &Nsig,NULL, NULL); 

printf(" csSDp=%.2E[cm~2] ==> Nsignal=°/ 0 .2E (TAB_1 column l)\n" ,csSDp[k] ,Nsig) ; 

} 

f_90=fluxFactorIC22( exLev,nu,nuB); 

printfC 907o°/ 0 exclusion limit for csSDp is °/ 0 .2E [cm~2] \n" , csSDp [k] *f _90) ; 
for(int i=l;i<NZ;i++) {nu[i]*=f_90; nuB[i]*=f_90;} 

IC22events(nu,nuB, 10,&Nsig,NULL,NULL); 

printfC 907J„ exclusion limit for number of signal events: °/ 0 .2E \n" ,Nsig/l. 2) ; 


The code for the calculation of an exclusion level in the framework of a given model 
reads 

double nu[NZ],nuB[NZ]; 
int forSun=l; 

WIMPSIM=1; 

neutrinoFlux(Maxwell,f orSun,nu,nuB); 

printf ("IceCube22 exclusion confidence level = °/„ .2E"/ 0 , /An" , 100*exLevIC22(nu,nu_bar ,NULL) ) ; 

For the default data hie mssmh.par provided in the MSSM directory of micrOMEGAs, 
we get a confidence level of 18%. 


C DarkSUSY 

DarkSUSY is generally used in the framework of some SUSY model. In order to compare 
the results obtained with micrOMEGAs in a model independent approach we use the fol¬ 
lowing trick to calculate the number of signal events for a specific channel and a DM 
- nucleon cross section. The DarkSUSY subroutine dsntICbounds provides the number 
of signal events expected in IceCube22 for given annihilation channels specified by the 
array wabr, for a DM mass wamwimp , and for DM nucleon cross sections with the pa¬ 
rameters wasigsdp for spin dependent and wasigsip for spin independent. We therefore 
redefine these parameters stored in common /wabranch/ before the call to the subroutine 
dsntICbounds. All cross sections should be specified in [cm 2 ]. The content of the ar¬ 
ray wabr is given in the hie DarkSUSY/src/wa/dswayieldone. f. In particular wabr(13) 
corresponds to the W + W~ annihilation channel and wabr(25) to the bb channel. 
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